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Abstract — In this work we seek to characterise the perfor- 
mance of spatial modulation (SM) and spatial multiplexing 
(SMX) with an experimental testbed. Two National Instruments 
(Nl)-PXIe devices are used for the system testing, one for the 
transmitter and one for the receiver. The digital signal processing 
that formats the information data in preparation for transmission 
is described along with the digital signal processing that recovers 
the information data. In addition, the hardware limitations of the 
system are also analysed. The average bit error ratio (ABER) 
of the system is validated through both theoretical analysis and 
simulation results for SM and SMX under line of sight (LoS) 
channel conditions. 

Index Terms — Spatial Modulation (SM), Spatial Multiplexing 
(SMX), Multiple-Input Multiple-Output (MIMO) systems, Ex- 
perimental Results, Wireless Testbed 



I. Introduction 

Multiple-input multiple-output (MIMO) systems offer a 
significant increase in spectral efficiency in comparison to sin- 
gle antenna systems [1,2]. An example is spatial modulation 
(SM), which increases the spectral efficiency of single antenna 
systems while avoiding inter-channel interference (ICI) [3]. 
This is attained as shown in Fig. 1, through the adoption 
of a new modulation and coding scheme, which foresees: 

i) the activation, at each time instance, of a single antenna 
that transmits a given data symbol (constellation symbol), and 

ii) the exploitation of the spatial position (index) of the active 
antenna as an additional dimension for data transmission (spa- 
tial symbol) [4]. Both the constellation symbol and the spatial 
symbol depend on the incoming data bits. An overall increase 
by the base-two logarithm of the number of transmit-antennas 
of the spectral efficiency is achieved. This limits the number 
of transmit antennas to be a power of two unless fractional 
bit encoded spatial modulation (FBE-SM) [5], or generalised 
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Fig. 1. The unique three dimensional constellation diagram for SM. The 
lower two bits, in the four bit word define the spatial-constellation point which 
identifies the active antenna. These are shown in parentheses. The remaining 
two bits determine the signal-constellation point that is to be transmitted. 

spatial modulation (GSM) [6] are used. In particular, in [6] it 
is shown that the number of spatial symbols does not need to 
be equal to the number of transmit antennas. For example, if 
GSM is used, the number of spatial symbols is equal to the 
number of unique channel signatures between the transmitter 
and receiver, where the unique channel signatures can be 
obtained by activating various combinations of the available 
transmit antennas. In this work, however, these unique channel 
signatures are assumed to be due to the activation of individual 
transmit antennas. 

Activating only one antenna at a time means that only one 
RF chain is needed, which significantly reduces the hardware 
complexity of the system [7]. Moreover, the most energy 
consuming part of a base station are the power amplifiers 
and the radio frequency (RF) chains associated with each 
transmitter [8], where the power requirements of a base station 
are shown to increase linearly with the number of RF chains 
added [9]. However, as only one RF chain is needed, SM 
offers a reduction in the energy consumption which scales 
linearly with the number of transmit antennas [10,11]. Fur- 
thermore, the computational complexity of SM-ML is equal 
to the complexity of single-input multiple-output (SIMO) 
systems [12], i.e. the complexity of SM-ML depends only on 
the spectral efficiency and the number of receive antennas, 
and does not depend on the number of transmit antennas. 
Moreover, in [13-15], the complexity of SM is further reduced 
by using the sphere decoder (SD). 
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Fig. 2. Block sequence of the main steps in the experiment, from the generation of the binary data to its recovery. 



Several papers that seeks to understand and hence improve 
the performance of SM in various scenarios are available 
in literature. In [16,17], the average bit error ratio (ABER) 
performance of SM is improved by introducing trellis coding 
on the transmitting antennas. The optimal detector is derived 
with and without channel state information at the receiver 
in [12, 18, 19]. The ABER performance is given when con- 
sidering channel estimation errors in [20-22]. The optimal 
power allocation for the case of two transmit antennas and one 
receive antenna system is given in closed form in [23], and 
the ABER performance of SM in correlated fading channels is 
considered in [24]. In [25,26] spectral efficiency and diversity 
gains are obtained by combining SM with space-time block 
codes (STBC-SM). Applying SM to relaying systems is 
also shown to result in significant signal-to-noise-ratio (SNR) 
gains when compared to non-cooperative decode and forward 
techniques [27]. In [28] the overall power performance of a 
base station (BS) employing SM is studied. More recently, a 
comprehensive analytical framework to compute the ABER 
of SM over generalized fading channels is introduced in [29] . 
Moreover, in [30] for the first time the performance of SM is 
analysed using real-world channel measurements. The latest 
research achievements and an outline of some relevant open 
research issues for SM are reviewed in [31]. All research thus 
far is strictly theoretical. 

In this paper, the ABER performance of SM is analysed in a 
practical testbed and compared with that for spatial multiplex- 
ing (SMX). In particular, the National Instruments (Nl)-PXIe- 
1075 chassis are used at the transmitter and receiver. The 
design of the testbed hardware and the software used are ex- 
plained in detail along with the transmission chain. The effects 
of the entire transmission chain on the system performance 
are examined. The basic elements of the transmission link are 
the transmit RF chain, the wireless channel, and the receive 
RF chain. In addition to the effects of the wireless channel 
on the phase and amplitude of the signal, the impact on the 
system performance of the power imbalances (Pis) in the 
transmitter and receiver RF chains is discussed. Furthermore, 
an analytical upper bound for the ABER performance of SM 
over non-line of sight (NLoS) channels with PI is derived, and 
compared to the experimental and computer simulation results. 
The experimental results validate the analytical bound as well 
as the attained computer simulations. Finally the performance 
of SM is compared with the theoretical and experimental 
results of SMX. 



This paper is organised as follows. The system set-up, 
equipment and digital signal processing are presented in 
Section II. The equipment constraints are then considered in 
Section EH, while the analytical modelling is discussed in Sec- 
tion IV. In addition, the computational complexity of the SM 
decoder algorithm is presented in Section V. The performance 
of SM is then characterised in the experimental and simulation 
environments in Section VI, where it is compared with the 
theoretical and experimental results of a SMX system. Lastly, 
the paper is summarised in Section VII. 

II. Testbed Set-up and System Model 

The testbed set-up and transmission chain can be separated 
into software and hardware parts, as shown in Fig. 2. The 
hardware consists of the Nl-PXIe chassis at the transmitter 
(PXIe-Tx) and the Nl-PXIe chassis at the receiver (PXIe-Rx). 
The software consists of the digital signal processing at the 
transmitter (DSP-Tx) and the digital signal processing at the 
receiver (DSP-Rx). 

The binary data to be broadcast is first processed by 
DSP-Tx, before being transmitted through the fading channel 
by the PXIe-Tx. The channel coefficient on the link between 
transmit antenna n t , and receive antenna r, is denoted by 
h( r nt y Note that the number of antennas at the transmitter 
and the receiver are denoted by Nt and N r , respectively. At 
the receiver, the PXIe-Rx records the RF signal and passes 
it through to the DSP-Rx for processing, where the original 
data stream is recovered. 

A. Testbed Hardware 

The NI-PXIe-1075 chassis are equipped with a 1.8 GHz 
Intel-i7 processor with 4 GB RAM and are shown in Fig. 3. 
The system has two transmit antennas and two receive anten- 
nas. Each antenna at the transmitter and receiver contains two 
quarter-wave dipoles, and one half-wave dipole placed in the 
middle. All three dipoles are vertically polarised. In addition, 
each antenna has a peak gain of 7 dBi in the azimuth plane, 
with an omnidirectional radiation pattern. 

1) Transmitter hardware (PXIe-Tx): The following NI- 
PXIe modules are used at the transmitter, 

. NI-PXIe-5450 16-Bit I/Q Signal Generator (SG-16bit), 

. NI-PXIe-5652 RF Signal Generator with a 500 kHz to 
6.6 GHz frequency range (SG-RF), 

• NI-PXIe-5611 intermediate frequency (IF) to carrier RF 
up-converter (up-converter). 



SERAFIMOVSKI el al: PRACTICAL IMPLEMENTATION OF SPATIAL MODULATION 



3 




(a) PXIe-Tx 



(b) PXIe-Rx 



Fig. 3. NI-PXIe-1075 chassis with the relevant on-board modules at the 
transmitter (PXIe-Tx), and at the receiver (PXIe-Rx). 



The PXIe-Tx has an operational frequency range of 
85 MHz to 6.6 GHz and can facilitate a bandwidth of 100 MHz 
at a maximum transmission power of 5 dBm. 

At the transmitter, the SG-16bit performs a linear mapping 
of the signed 16-bit range to the output power and polarisation, 
i.e., the peak voltage amplitude is assigned to any value in 
the transmission vector equal to 2 15 with a linear scale of 
the voltage amplitude down to zero. The output from SG- 
16bit then goes to SG-RF, which is connected to the up- 
converter. The up-converter outputs the analogue waveform 
corresponding to the data resulting from DSP-Tx at a carrier 
frequency of 2.3 GHz. This completes a single RF chain. 
The transmission of the RF signal by the up-converters is 
synchronised by using a 10 MHz reference signal. 

2) Receiver hardware (PXIe-Rx): The following Nl-PXIe 
modules are used at the receiver, 

. NI-PXIe-5652 an on-board reference clock (SG-RF), 

• NI-PXIe-5622 16-Bit Digitiser which records data sam- 
ples in an 116 format (16-Bit Digitiser), 

• NI-PXIe-5601 RF down-converter (down-converter). 
The PXIe-Rx can operate in a frequency range of 10 MHz 

to 6.6 GHz and can facilitate an operational bandwidth of 
50 MHz. For more details about the specifications of each 
module the reader is kindly referred to [32,33]. 

At the receiver, each antenna is associated with a complete 
RF chain. For each antenna, the down-converter is used to 
detect the analogue RF signal from its dedicated antenna. The 
signal is then sent to the dedicated 16-Bit Digitiser. The 16- 
Bit Digitiser applies a bandpass filter with a real flat bandwidth 
equal to Bf = (0.4 x f s ), where f s is the sampling rate 
[32]. The sampling rate in the experiment is 10 Ms/s which 
results in a real flat bandwidth of 4 MHz. This may result 
in frequency-selective fading. Nonetheless, equalisation is not 
required for the detection of SM or SMX signals in this 
experiment because: i) there are no multi-tap delays in the 
experimental setup due to very small distance between the 
transmit and receiver antennas; and ii) maximum likelihood 
(ML) detection is used to decode the receiver signal for both 
SM and SMX. The use of ML detection is applied to the 
complete SM symbol, i.e., the spatial symbol and the signal 
symbol are decoded jointly. Finally, after synchronisation of 
the 16-Bit Digitiser with the on-board reference clock of the 
SG-RF, the 16-Bit Digitiser writes the received binary files. 
The simultaneous recording of the two signals coming from 
Txl and Tx2 is facilitated by utilizing multiple processing 
cores and multiple Nl-PXIe modules. The recorded files are 
then processed according to DSP-Rx in Fig. 4. 



DSP-Tx 



Sync and SNR 



Tuning 
Signal 
Power 



Up Sampling and 
Filtering 



Pilot and Zero 
Padding 



SM Modulation 



Framing 



Binary file for Tx1 

^\ I Binary file from Rx1 




DSP-Rx 

- Synchronisation 

SNR Calculation 

Extract Frames 

Down Sampling and 
Filtering 



Frequency Offset 
Correction 

T 



Channel 
Estimation 



Binary Data 



SM De-Modulation 



Binary Data 



Fig. 4. A step-by-step layout of the binary data encoder (DSP-Tx) and 
decoder (DSP-Rx) processes. 



B. Testbed Software 

Matlab was used to facilitate the digital signal processing re- 
quired at the transmitter, DSP-Tx, and at the receiver,DSP-Rx. 
DSP-Tx processes the information data and generates binary 
files that can be transmitted by PXIe-Tx. DSP-Rx process the 
received data from PXIe-Rx and recover the original informa- 
tion data stream. Fig. 4 outlines the processing algorithms at 
DSP-Tx and DSP-Rx. 

1) DSP-Tx: The DSP-Tx process takes the incoming 
binary information data and performs the following, 

1.1 Framing: The incoming data is split into frames consist- 
ing of 100 symbols per frame. 

1.2 Modulation: The data in each frame is modulated using 
SM or SMX: 

• SM: The bit stream is divided into blocks containing 
log 2 (Nt M) bits each, where M is the signal constella- 
tion size. The following mapping rule is then used [4]: 

a) The first log 2 (Nt) bits determine which transmit 
antenna is active, i.e., they determine the spatial 
constellation point of SM. In this paper, the trans- 
mit antenna broadcasting is denoted by n t with 
n t e{l,2,...,N t }. 

b) The second log 2 (M) bits are used to choose a 
symbol in the signal-constellation diagram. Without 
loss of generality, quadrature amplitude modulation 
(QAM) is considered. The actual complex symbol 
emitted by the transmit antenna n t is denoted by 
s t , with s t € {s 1 ,s 2 ,...,s M }- 

By following the above steps, the JV( x 1 dimensional 
transmit vector is: 

x n t ,s t = [Oix(n t -l), S t , Olx(Nt-nt)] > (1) 

where [■] T denotes the transpose operation, and pxq 
is a p x q matrix with all-zero entries. Equation (1) 
is a representation of the transmission vector for SM. 
Since SM activates only one transmit antenna at any 
transmission instance, it means that only one transmit 
antenna can broadcast a symbol while all others remain 
silent. To this extent, the transmit vector is composed 
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of all zeros, except for the single symbol, St, which is 
broadcast from antenna n t . In this manner, SM avoids 
ICI and allows single-stream ML decoding. In addition, 
SM is energy efficient since only a single RF chain is 
active, while still providing a multiplexing gain [10]. 
• SMX: In this case, the bit stream is divided into blocks 
of iV t log 2 (M) bits, then, according to [34]: 

a) Each log 2 (M) bits are separately modulated using 
M-QAM modulation. 

b) The modulated symbols are then transmitted simul- 
taneously from the N t transmit antennas. 

1 .3 Pilot and Zero Padding: The least squares (LS) channel 
estimation algorithm with local orthogonal pilot sequences 
is used to estimate the channel [35]. Two pilot signals are 
added for each frame, one at the start of the frame, and one 
at the end. Each pilot signal contains ten pilot sequences, 
where the orthogonal pilot sequence for the n t -th transmit 
antenna is defined as, 



(2) 



where 0, H (£) is the ^-th element of the pilot sequence 
O nt transmitted from antenna n t , j = \/— 1 is the 
imaginary unit and Ne is the cardinality of the pilot 
sequence. In this work, the length of each pilot sequence 
is Nq = 10. To avoid inter-frame interference (IFI), an all 
zero sequence of 50 zero valued symbols is added to both 
the start and the end of the frame. Furthermore, a sequence 
of constant valued symbols is added to enable frequency 
offset (FO) estimation at the receiver. The length of the 
FO estimation sequence is 1000 symbols. 

1.4 Up Sampling and Filtering: Up-sampling and matched 
filtering (pulse shaping) are used to maximise the SNR 
and reduce inter-symbol interference (ISI) [36]. Each 
frame is up-sampled with an up-sampling ratio of 4, 
and then passed through a root raised cosine (RRC)-finite 
impulse response (FIR) filter with 40 taps and a roll-off 
factor of 0.75. The large roll-off factor is necessary to 
ensure that the power is focused in a short time instance 
to ensure that only a single RF chain is active when using 
SM. 

1.5 Tuning Signal Power: The SNR is varied by changing 
the power of the transmitted signal to obtain the ABER. 
This is done by multiplying each transmission vector with 
a "Tuning Signal Power" factor to obtain the desired 
transmit power. In particular, by changing the amplitude 
of the "Data section" in the transmission vector by using 
the "Tuning Signal Power" factor. 

1.6 Synchronisation and SNR: Several preamble- 
autocorrelation based methods for frame synchronisation 
were tested [37-39]. However, despite the introduction 
of an interpolation filter at the receiver and due to 
the channel attenuations, the estimated start of the 
signal was typically in error by one or two samples. 
This meant that sample synchronisation could not be 
achieved consistently, resulting in off -by-one errors. The 
investigation of synchronisation techniques is outside the 



scope of this work, but in order to avoid synchronisation 
via a cable, as is often done in similar experimental 
systems, the peak detection technique has been applied 
which resulted in the desired outcome. We recognise 
that this technique is suboptimal as it results in power 
amplifier saturation and potential signal distortions. 
Nonetheless, a sequence of 20 symbols with maximum 
power, separated by 50 zero valued symbols between 
each, are added to the start of the transmitted signal. 
The large power difference between the maximum power 
peaks and the power of the "Data section" symbols is 
reasonable since the instantaneous channel power may 
fluctuate by as much as 20 dB due to fast fading [40, 
41]. The power difference between the synchronisation 
section and the remaining sections is set to be larger 
than the maximum channel variation. In this manner, a 
successful peak detection is guaranteed. If this is not the 
case, no peak may be detected at the receiver and all 
further decoding would be erroneous. 
To facilitate SNR calculations at the receiver, two se- 
quences of power and no power are added after the 
synchronisation pulses of the transmitted signal, indicated 
by "SNR section" in Fig. 5. Each sequence contains 
5 blocks of 50000 symbols and 50000 zeros. The first 
sequence is transmitted from the first antenna while the 
second antenna is off. The second sequence is transmitted 
from the second antenna while the first antenna is off. 
After the DSP-Tx process completes, the transmit vector 
symbols are converted to 116 format and are recorded to a 
binary file. This binary file is then broadcast by PXIe-Tx. 

Fig. 5 is an absolute value representation of the processed 
incoming data that is passed to the first transmit antenna (Txl) 
and Fig. 6 shows the absolute value representation of each 
frame. Note that the "Data section" is a series of concate- 
nated frames. In Fig. 6, it can be seen that each frame 
contains 26100 samples. Therefore, the period of each frame 
is Tprame = 26100// s = 2.6 ms. This is much less than the 
coherence time of the channel given that, typically, the co- 
herence time for a stationary indoor environment is approxi- 
mately 7ms [41, and references therein]. Hence, the channel 
estimation at the receiver is valid for the frame duration. 

2) DSP-Rx: The data received by PXIe-Rx is processed 
by DSP-Rx to recover the original data stream. To accomplish 
this, the following steps are required: 

2.1 Synchronisation: This is achieved by searching for the 
peaks with a value above a certain threshold in the 
received signal. The threshold is set as 70% of the highest 
value in the received vector. This threshold level accounts 
for the natural voltage variations in the system, i.e., the 
difference between peak voltage and root-mean-square 
voltage. If the number of peaks found is less than 20, then 
the received vector is discarded from further calculations. 

2.2 SNR Calculation: The SNR is defined as, 



E 



SNR 



IHxII 



(3) 



where H is the N r x Nt channel matrix, x is the N t x 1 
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Fig. 5. This is the absolute value representation of the transmission vector 
being sent to Txl. The synchronisation, SNR estimation and Data sections 
are shown. The value of the peak must equal 2 lj since the 16bit-Digitiser 
operates using an 116 format before tuning the signal power of the data. The 
highest value in the SNR section is the same as the highest value in the 
information data section, in this example a value of 2896. The peak value is 
2 15 . There is approximately a 21.1 dB difference between the peak power in 
the synchronisation section and the peak power in the SNR estimation and 
data sections. This is apparent when looking at the two data points shown in 
the figure. 
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Fig. 6. This is the absolute value representation of a single frame from the 
vector being transmitted by Txl in the 116 data format, which is a signed 
15 bit representation of an integer number. 



transmitted vector, E [•] is the expectation operator, and 
|| • || F is the Forbenius norm. 

Assuming that the noise at the receiver is additive white 
Gaussian noise (AWGN), the received signal for the 
duration of the SNR sequence can be written as follows: 



y = h„ t s t + n 



(4) 



where y is the N r x 1 received vector, h„ t is n t column of 
the channel matrix H, n is the iV r x 1 AWGN vector with 
cr^ variance and /i„ mean, and s t is the transmitted symbol 
from the n t antenna. As mentioned in Section II-B1, only 
a single transmit antennas is active when broadcasting 
the SNR sequence and s t is either equal to the maximum 
value in the "Data section" x lmx or zero, as shown in 



Fig. 5. Hence, the received signal in (4) can be re-written 

as, 

hn t ^max ~t~ 1; &t — -^max 

n, s t = 



(5) 



Proceeding from (5), 
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IHxIl 



E 



= E 



n 



E[||n||f 



(6) 
(7) 



where [•]■" is the Hermitian operation. As discussed in 
Section II-B1, each SNR sequence contains 50000 sym- 
bols and 50000 zero valued symbols. Since the noise in 
the system represents an ergodic process, the ensemble 
average in (6) can be replaced with a time average, 



E 



IHxIl 



E (WIp- Wp-2yfni 



.-,()! MID 



E II* 



»IIF 



.,(11101) 



E IN 



i=l 



(8) 
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where y, and are the z— th received vector. To get a 
more accurate estimation, the SNR is calculated for the 5 
transmitted SNR sequences received at both antennas and 
then averaged again over those measurements. 

2.3 Extract Frames: After finding the start of the transmis- 
sion and calculating the SNR, DSP-Rx performs a serial 
to parallel conversion to separate the received frames. 

2.4 Down Sampling and Filtering: To complete the matched 
filter described in Section II-B1, each frame is down- 
sampled by a factor of 4 and passed through an RRC-FIR 
filter. 

2.5 Frequency Offset (FO) Correction: The DSP-Rx esti- 
mates the FO for each frame by, 

Zxinrin — ^-Xt, 

A f = — — (10) 

3 2tt x 1000 

where Zxiooo, Zxi are the angles of the first and the last 
sample of the FO sequence transmitted by the DSP-Tx 
where the FO sequence has exactly 1000 symbols. These 
angle values are obtained by correcting the radian phase 
angles in a vector by adding multiples of ±27r as required. 
This enables a better estimate of the phase offset. As- 
suming a linear phase rotation, the frequency offset can 
be estimated using (10). The FO for each frame is then 
corrected by, 



y, 



Vi x e 



— j"27T A fi 



(ii) 



where j/j , y; are the i-lh element of the corrected and the 
uncorrected received frame, respectively. 
2.6 Channel Estimation: The channel estimation is done by 
using the LS channel estimation algorithm proposed in 
[35], where for each frame the channel is estimated by, 

1 



H 



LS 



N e 



H H r 



(12) 



where H r is the received pilot sequence. To enable a more 
accurate evaluation of the system, the channel is estimated 
and averaged over 10 pilot sequences. Furthermore, two 
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Fig. 8. Physical experimental layout: A pair of receive and a pair of transmit antennas are set 2.2 m apart from each other with a direct line of sight. Each 
pair of antennas is 1.5 m from the ground and there is a 10 cm spacing between the antennas in either pair corresponding to 0.77 times the wavelength at 
2.3 GHz. All antennas are omnidirectional. 



channels are estimated per frame, the first channel esti- 
mate is used for the first half of the data symbols in the 
frame, and the second is used for the second half of the 
data symbols in the frame. 
2.7 Demodulation: The ML optimum receiver for MIMO 
systems is used, which can be written as, 



-(ML) 



arg mm 

x£ Q 



in{||y-Hx||*) (13) r 



where Q contains every possible (N t x 1) transmit vector, 
and : denotes the estimated transmission vector. However, 
since only one transmit antenna is active at a time for a 
SM system, the optimal receiver in (13) can be simplified 
to, 



N r 



(ML) .(ML) 



arg mm 

n t e{i,2,...,JV t } 

s€{si,s 2 ,...,s M } 



l (r,n t )-- 



(14) 

where y r is the r-th entries of y. 
Finally, the recovered binary data along with the estimated 
SNR are used to obtain the ABER performance of both SM 
and SMX. 
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Fig. 9. CDFs for each of the fast fading coefficients, h(r,nt)> °f me rour 
channels in the experiment. Each is defined by a Rician distribution with a 
unique i\"-factor. The markers denote the measurement points while the lines 
denote the best fit approximation. Note that the wireless channel mean values 
fall in the range of 1.3 mV to 3.6 mV. 



C. Propagation Environment ( Channel) 

The physical layout of the experimental set-up is shown in 
Fig. 7 and the relative antenna spacing is provided in Fig. 8. 
In particular, the two transmit and two receive antennas are 
identical and are placed directly across from each other. As 
such, the channel between the transmitter and receiver has a 
strong line of sight (LoS) component. Therefore, the channel 
is assumed to be a Rician fading channel with a large K- 
factor due to the short distance between the transmit and 




Fig. 7. Experimental setup in the laboratory. 



receive antennas, where K is the ratio of the coherent power 
component, usually the LoS, to the non-coherent power com- 
ponents, usually NLoS. The omnidirectional transmit antennas 
broadcast on a frequency of 2.3 GHz at 10 Ms/s. 

Channel measurements were collected to verify that the 
channel environment followed a Rician distribution. To 
achieve this, the transmitter broadcasts pulses at 10 Ms/s on 
a carrier frequency of 2.3 GHz at 4 dBm peak power. Each 
pulse includes a frequency offset estimation section and a total 
of 10 5 pulse samples are collected. A best fit approximation is 
then calculated for the collected data. In particular, a maximum 
likelihood estimation is fitted to the collected data. A Chi- 
squared goodness-of-fit test is then performed to ascertain 
that the distribution resulting from the maximum likelihood 
estimation fits at least 95% of the data. The empirical CDF 
for each link is presented in Fig. 9. The results show that 
the channel does follow a Rician distribution with a K factor 
that ranges between 31 — 38 dB. The different i\~-factors on 
the links between the transmit and receive antennas may be 
explained by the room geometry, the antenna positioning and 
the overall propagation environment. However, note that each 
of the CDFs has a different mean, which will be discussed in 
the next section. 
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Fig. 10. CDFs for each of the fast fading coefficients, h( r ,nt)> °^ tne ' our cnanne l s in me experiment. Each is defined by a Rician distribution with a 
unique iiT-factor. The markers denote the measurement points while the lines denote the best fit approximation. Despite using a coaxial cable with a 10 dB 
attenuation to connect the RF chains, each channel exhibits a unique mean. 



III. Equipment Constraints 

Fig. 8 shows the physical layout of the experiment. Note 
that the 10 cm inter-antenna separation used here is suffi- 
cient to guarantee very low, if any, spatial correlation when 
broadcasting at 2.3 GHz with a 2.2 m separation between the 
transmitter and receiver [42]. 

The physical environment through which the signal passes, 
starting from the SG-RF at the transmitter, until it reaches 
the 16-Bit Digitisers at the receiver, suffers from connector 
losses, differences in the RF chains, different phase responses, 
attenuations and similar. To study and model the effects of the 
hardware imperfections on the signal power: 

• An RF coaxial cable with a 10 dB attenuation is con- 
nected between each transmit and receive antenna. 

• A pulse is transmitted at 10 Ms/s on a carrier frequency 
of 2.3 GHz at —10 dBm peak power. Each pulse includes 
a frequency offset estimation section and a total of 10 5 
pulse samples were collected. 

• The CDF for each of the fading coefficients is calculated 
and is shown in Fig. 10. 

In an ideal environment, the means of the CDFs in Fig. 10 
should be equal. However, imperfections in the hardware result 
in different means for each transmit to receiver antenna pair, as 
can be seen in Fig. 10. The differences between the channels 
can be modelled as a PI between the various link pairs in 
the channel matrix H. Therefore, the channel coefficients are 
redefined as, 



where oef r<n A is the channel attenuation coefficient from 
receive antenna r to transmit antenna n t . 

To locate the source of the discrepancy between the different 
channel attenuations, i.e., determine if the NI modules or the 
NI chassis is the source, the RF chains at the receiver were 
swapped around and the channels were estimated in config- 
uration (I) and configuration (II). To clarify, configuration 
(I) represents the default modular set-up of the testbed while 
configuration (II) refers to swapping the front-end modules 
around the transmit chassis. Fig. 10(a) shows the channel CDF 
for each transmit to receive antenna pair in configuration (I) 
while Fig. 10(b) shows the channel CDF for each transmit 
to receive antenna pair in configuration (II). By considering 
the means of the CDFs in Figs. 10(a) and 10(b) and taking 
/i( L1 j as a base, the various channel attenuations that result 
when the receiver is in configuration I or in configuration II 
are given in (16) and (17) respectively. Comparing Fig. 10(a) 
and Fig. 10(b), as well as the attenuations in (16) to those in 
(17), shows that they are very similar. Indeed, swapping of 
the RF chains has a minimal impact on the estimated mean 
of each channel attenuation. Thus, it can be assumed that 
the NI modules that compose the receive RF chains are the 
source of the hardware imperfections, and consequently lead 
to the differences in the means of the estimated CDFs. To 
account the hardware imperfections, the channel attenuation 
coefficients in (16) and (17) are taken into consideration in the 
derivation of the analytical model in Section IV. The accuracy 
of the derived analytical bound using the channel attenuation 
coefficients in (16) and (17) is demonstrated in Section VI 
where it is compared with empirical results. 



o (M) =0dB, a (2) i) = 0.25 dB, a {h2) = 0.88dB, a (2)2 ) = l.ldB, 
a ( i )1) =0dB, a (2 ,i) = 0.29dB, a (1 , 2 ) = 1.13dB, a (2 , 2 ) = 1.17 dB. 



(16) 
(17) 
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IV. Analytical Modeling 

An analytical model for the ABER performance of the 
experimental system is developed by considering the system 
model presented in Section II and the system constraints in 
Section III. The performance of SM and SMX over a single 
link in a noise-limited scenario is characterised by 

ABER <iVV ^il^E„ jPr), (18) 

- 2"i jL^, m n tenor] 

Xf x 

where N (x t , x) is the number of bits in error between 
the transmitted vector x t and x, Eh{-} is the expectation 
across the channel H, and Pr is the conditional pairwise 

error 

error probability (PEP) of deciding on x given that x t is 
transmitted [43], 
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Fig. 11. Relative receiver complexity reduction of the SM-ML receiver 
versus the SMX-ML receiver. 

. SM-ML: The computational complexity of the SM-ML 
receiver outlined in (14) is equal to, 



Pr 



Pr( lly-HxtH' > ||y-Hx||| 



H 



Csm-ml — 8N r 2 r 



(22) 



Q V7ex||H(x t -x) 



(19) 



tH^- is half of the SNR between the transmitter 

21\q 



where 7 c; 

and receiver, and Q(oj) — /J° exp (— 4^ dt is the Q- 
function. As Fig. 8 indicates, the transmit and receive antennas 
in the experiment experience a very strong LoS environment. 
Accordingly, the channel between each transmit to receive 
antenna pair is characterised by Rician fading. A generic 
Rician channel is defined as 



K 



1 + K 



1 



1 + K 



l (r,n t ) 



(20) 



where hu n A ~ CAf(0, 1) is a complex normal, circular sym- 
metric random variable with zero mean and unit variance. 
n t € {1, 2} is the index of the transmit antenna and re {1,2} 
is the index of the receive antenna. 

To account for the hardware imperfections that result from 
the power imbalances, the fast fading channel coefficients 
are redefined according to (15), (16), and (20). Section VI 
validates the derived analytical bound by comparing it to 
experimental and simulation results. 

V. Complexity Analysis 

The computational complexity of the ML detector for SM 
(SM-ML) is compared to that of the ML detector for SMX 
(SMX-ML). The complexity is computed as the number of 
real multiplicative operations (x,-=-) needed by each algo- 
rithm. The detailed derivation of each expression is considered 
in [13] and references therein. 

* SMX-ML: The computational complexity of the SMX- 
ML receiver outlined in (13) is equal to, 



4(N t + l)N r 2 m , 



(21) 



where m is the spectral efficiency of the system. Note 
that ^|y — Hx| 2 ^j in (13) requires (Nt + 1) complex 
multiplications. 



where the ML detector searches through the entire trans- 
mit and receive search spaces. Note that evaluating the 
Euclidean distance (\y r — h^ rnt yst\ 2 ^j requires 2 com- 
plex multiplications, where each complex multiplication 
requires 4 real multiplications. 

Considering (21) and (22), for the same spectral efficiency, 
the reduction in complexity of the SM-ML receiver relative 
to that of the SMX-ML receiver is given by, 



C rd = 100X1 



N t + 1 



(23) 



On the one hand, as can be seen in (22), the complexity of 
the SM receiver does not depend on the number of transmit 
antennas, and it is equal to the complexity of SIMO systems. 
On the other hand, the complexity of SMX increases linearly 
with the number of transmit antennas. Therefore, as the 
number of transmit antennas increases, the relative gain of 
the SM receiver increases. This can be seen in Fig. 1 1 where 
the relative complexity for Nt G {2, 4, 8, ... , 128} is shown 
for both systems. In fact, Fig. 1 1 shows that for Nt — 4, SM 
offers a 60% reduction in complexity, while a 98% reduction 
in complexity can be seen for Nt — 128. The theoretical, 
simulation and empirical results for SM and SMX are now 
discussed. 

VI. Experimental Results and Numerical Analysis 

A. Measurement Campaign 

A stream of 10 5 information bits is sent per transmission 
to obtain the experimental results. Two transmit antennas are 
available and binary phase shift keying (BPSK) is used for 
the signal constellation. As mentioned in Section II-A2, the 
real flat bandwidth is 4 MHz. The information data is put in 
50 frames with 2000 bit each, as shown in the "Data section" 
of Fig. 5. The channel is estimated at the beginning and the 
end of every frame, resulting in 100 channel estimations per 
transmission vector. The experiment is repeated 1000 times for 
every SNR point. In addition, analytical and simulation ABER 
curves are shown for SM in a Rician environment with and 
without the Pis given in (16). 
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B. Results 

The simulation, analytical and experimental results for the 
ABER performance of SM in a LoS channel are illustrated in 
Fig. 12. In particular, the experimental results approximate the 
performance of the simulation results with Pis and both the 
simulation, and experimental results, are closely approximated 
by the derived upper bound at a low ABER. 

This result serves to validate theoretical work done in the 
field where the presented SNR along the x-axis is equivalent 
to the SNR on The large error between the experi- 

mental, simulation and analytical curves at high ABER can 
be attributed to a number of factors including incorrect fre- 
quency offset estimation, timing recovery errors, synchronisa- 
tion problems, poor channel estimation and decoding. Notably, 
incorrect frequency offset estimation can result in a systematic 
error contributing significantly to the 30% error seen at low 
SNRs in the figure. As the SNR increases, however, frequency 
offset estimation, timing recovery and channel estimation 
improve, leading to a lower ABER as shown in Fig. 12. Dif- 
ferences between the measured and simulated ABER curves 
can be attributed to channel imperfections such as channel 
correlations, mutual coupling and interference signals from 
the surrounding environment. Quantifying these imperfections 
is deemed important and requires channel modelling and 
interference measurement. However, addressing these effects 
is beyond the scope of this work and will be subject of future 
works. 

SM performs best in a rich scattering environment where the 
channel between each transmit and receive antenna is unique. 
In particular, the larger the Euclidean distance between two 
received vectors is, the better the ABER performance of SM 
becomes. Conversely, the more similar the channels are, the 
worse the ABER of SM is. However, the channel uniqueness 
can be the result of the scattering environment or Pis caused 
by hardware tolerances. The analytical and simulation results 
presented in Fig. 12 show the poor performance of SM 
in a Rician environment with no PI between the various 
transmitter to receiver links. Fig. 12 also shows the analytical 
and simulation ABER for SM when PI are introduced. Indeed, 
the ABER of SM improves significantly when these Pis are 
introduced as each channel becomes more separable. This 
increases the Euclidean distance and improves performance. 

If the channels between each transmit antenna to each 
receive antenna are similar, then the ABER performance of 
SM degrades. This is seen when looking at the SM system 
without Pis, illustrated by the dashed green line with trian- 
gular markers in Fig. 12. In fact, the ABER of SM can be 
approximated by separating the error that originates from the 
estimation of the spatial constellation symbol and the error 
that originates from the estimation of the signal constellation 
symbol [44]. Therefore, depending on the environment, the 
main contributor to the overall ABER of a SM system will be 
the erroneous detection of the spatial or signal constellation. 

When Pis are introduced, the Euclidean distance between 
the channel signatures increases. This decreases the error 
contribution of the spatial component of SM. Hence, when 
the SNR is sufficiently high to have near perfect channel 
estimation, the error of the system is bound by the error from 
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Fig. 12. ABER for SM in an experimental set-up with 2 transmit antennas, 
2 receive antennas and a spectral efficiency of 2bits/s/Hz. The SNR is set 
as measured on fen i) with <xn i) = OdB. The solid black line with square 
markers denotes the experimental results. The green diamond markers denote 
simulation results with no power imbalance (PI) between the links while the 
green dashed line is the analytical prediction. The remaining curves denote 
the simulation (Sim) and analytical (Ana) results. 

the signal component of SM. This separation can only be 
shown when iterative detection is used, which is proven to 
be sub-optimal [12]. In addition, work in [45] shows that the 
error when only the spatial constellation of SM is used for data 
transmission gets worse for an increasing K factor in a Rician 
environment. This is the opposite to conventional modulation 
techniques since a larger K factor for SM means a smaller 
Euclidean distance between the spatial constellation points 
which results in an increased ABER performance. Indeed, 
it is the Euclidean distance between the different channels 
that determines the error in the spatial constellation detection. 
However, since ML-optimal detection is used at the receiver, 
separating the error from the spatial and signal symbols is 
strictly not permitted. Please note that the Pis between the 
links are always obtained relative to the channel with the 
greatest attenuation, i.e., the values of the PI factors in (16) 
and (17) are always positive. 

Furthermore, power imbalances between the transmitting 
antennas are shown to offer improved performance in terms of 
the ABER when only the spatial constellation of SM is used, 
i.e., when space shift keying (SSK) is the underlying modu- 
lation technique. In particular, an optimised power allocation 
for a various number of transmit antennas is addressed in [23], 
where the authors show that there is optimal power allocation 
between the transmitting antennas which can serve to increases 
the Euclidean distance between the channel signatures and 
improve the ABER performance of SM. Indeed, SM has 
also been successfully applied to an AWGN optical wireless 
channel where it is shown that Pis greatly improve the ABER 
performance [46]. 
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Fig. 13. ABER for SMX in an experimental set-up with 2 transmit antennas, 
2 receive antennas and a spectral efficiency of 2bits/s/Hz. The SNR is set 
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markers denotes the experimental results. The green diamond markers denote 
simulation results with no power imbalance (PI) between the links while the 
green dashed line is the analytical prediction. The remaining curves denote 
the simulation (Sim) and analytical (Ana) results. 



The simulation, analytical and experimental results for the 
ABER performance of SMX in a LoS channel are illustrated 
in Fig. 13. In particular, the experimental results closely follow 
the performance of the simulation results with Pis and both the 
simulation, and experimental results, are closely approximated 
by the derived upper bound at low ABER when the hardware 
imperfections are taken into account. This result serves to 
validate theoretical work done in the field. The results in 
Fig. 13 demonstrate that the SMX system, like the SM system, 
also benefits from the Pis in the hardware. The SMX system 
exhibits approximately a 3 dB coding gain when compared 
to SM at an ABER of 10~ 4 . This coding gain can also be 
seen at an ABER of 10~ 3 in Fig. 14, where the simulation 
and analytical results for the ABER performance of SM and 
SMX are shown when there are no Pis between the links. 

The coding gain of SMX relative to SM is expected when 
there are few transmit antennas. The Euclidean distance be- 
tween the transmit vectors, and therefore the variance in (19), 
in SMX is larger than in SM. However, the aim of this 
paper is to show that empirical results validate the simulation 
and analytical work done in the field, which can be seen in 
both Fig. 12 and Fig. 13. Unfortunately, due to the limited 
number of transmitter and receiver RF chains available, there 
are no experimental results for systems with a larger number 
of transmit or receive antennas where SM is shown to perform 
better than SMX. These empirical results will be the focus of 
future research. Nonetheless, the accuracy of the theoretical 
and simulation results of SMX and SM with a large number 
of transmit and receive antennas can be extrapolated from the 
presented results. 



Fig. 14. ABER for SM and SMX in a Rician fading channel where 
K = 33 dB with 2 transmit antennas, 2 receive antennas, a spectral efficiency 
of 2bits/s/Hz and no Pis between the channels. The simulation (Sim) are 
denoted by markers while the analytical (Ana) results are denoted by the 
lines. 
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Fig. 15. Simulation results for the ABER for SM and SMX in a Rayleigh 
fading environment with a spectral efficiency of 8 bits /s/Hz and no Pis 
between the channels. 



Fig. 15 compares the ABER between SM (solid lines) and 
SMX (dashed lines) in a system with a large number of 
transmit antennas. Each system operates in a Rayleigh fading 
environment with a spectral efficiency of 8 bits/s/Hz and four 
receive antennas. The results demonstrate the coding gains 
available to a SM system as compared to SMX when a large 
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number of transmit antennas are available. In particular, SM 
with 7V t = 64 offers a coding gain of up to 4 dB with 
respect to SMX with N t — 8 and a coding gain of 6 dB 
with respect to SMX with N t — 4. These performance gains 
stem from the greater Euclidean distance between the transmit 
vectors for SM. It is important to note that although SM is 
simulated as having 64 transmit antennas available, it requires 
only a single RF chains, while SMX requires 8 RF chains for 
the 8 transmit antennas. Furthermore, to achieve the ABER 
performance illustrated in Fig. 15, SM requires 64 unique 
channels. In this work, a unique channel is assumed to be 
available only with the addition of a single transmit antenna. 
However, work in [6, 25] and others, look at creating multiple 
channel signatures without the need for a large number of 
physical transmit antennas while maintaining a similar ABER 
performance to the traditional SM scheme. 

This work demonstrates that the hardware tolerances of 
practical communication systems are beneficial for the ABER 
performance of both SM and SMX. This behaviour along with 
the requirement for a single RF chain, make SM a viable 
candidate for future wireless networks. 

VII. Summary and Conclusion 

In this work, the ABER performance of SM and SMX has 
been validated experimentally for the first time. In particular, 
the encoding and decoding processes were presented. The 
experimental testbed, equipment and channel conditions were 
then described in detail and the ABER of SM and SMX were 
obtained in a practical testbed environment. In addition, the 
experimental results were compared to both simulation and 
analytical approaches. As a result, it has been shown that 
a Rician channel with different channel attenuations closely 
described the behaviour of SM and SMX in the physical envi- 
ronment. Furthermore, it was demonstrated that the different 
channel attenuations resulted from various hardware imper- 
fections at the transmitter and receiver RF chains. In fact, the 
induced power imbalances resulted in significant coding gains 
for the practical systems relative to the theoretical predictions 
without such power imbalances. To this extent, SM and SMX 
performed as expected relative to the theoretical work when 
the power imbalances were introduced in the analytical model. 
This result validated the SM principle. The performance gains 
exhibited by SM in the practical implementation make SM a 
viable candidate for future wireless networks and particularly 
for systems with a large number of transmit antennas available. 

It is worth noting that the presented work may be extended 
in a number of different ways that would broaden its applica- 
bility. Empirical results that demonstrate the performance of 
SM and SMX with a large number of transmit and receive an- 
tennas remain to be obtained. In light of the above results, the 
ABER performance of SM and SMX is expected to follow the 
theoretical models, but these results are essential to validate 
the ABER performance for both SM and SMX systems. In ad- 
dition, channel imperfections such as channel correlations and 
mutual antenna coupling along with the impact of interfering 
signals from neighbouring transmitters on the same frequency, 
should be analysed. Furthermore, obtaining empirical results 
for the capacity and energy efficiency of SM are of great 
interest for future research, particularly since SM is proj ected 



to have large energy efficiency gains when compared to other 
traditional MIMO schemes since it requires only a single RF 
chain. As a consequence, the quiescent power and circuit 
power can be kept at low levels. Acquiring the hardware which 
would enable the accurate measurement of these aspects is key. 
Lastly, the implementation of the SM detection algorithm on 
a DSP or an FPGA brings with it a number of optimisation 
challenges such as the use of multi-threading, pipeling, fixed 
point computations and others. The deployment of SM on an 
FPGA or a DSP has yet to be demonstrated. 

It has been shown that SM is a simple, low cost, MIMO 
technique, which has now demonstrated excellent performance 
in a LoS wireless channel. Therefore, this work shows that SM 
is a promising practical approach to obtaining the enhanced 
performance of spatial multiplexing without introducing high 
processor complexity and high power consumption that would 
occur when using other spatial multiplexing approaches. The 
aim now is to investigate the performance of SM in a range of 
experimental channel conditions and further study its potential. 
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